Model of tunneling transistors based on graphene on SiC 
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Recent experiments shown that graphene epitaxially grown on Silicon Carbide (SiC) can exhibit 
a energy gap of 0.26 eV, making it a promising material for electronics. With an accurate model, 
we explore the design parameter space for a fully ballistic graphene-on-SiC Tunnel Field-Effect 
Transistors (TFETs), and assess the DC and high frequency figures of merit. The steep subthreshold 
behavior can enable /on/^off ratios exceeding 10 4 even with a low supply voltage of 0.15 V, for 
devices with gatelength down to 30 nm. Intrinsic transistor delays smaller than 1 ps are obtained. 
These factors make the device an interesting candidate for low-power nanoelectronics beyond CMOS. 



Carbon-based nanoelectronics is characterized by rel- 
atively small bandgaps, and symmetric dispersion rela- 
tions for electron and holes, associated to a relatively 
small effective mass. This is the ideal situation for 
TFETs, whose carbon implementation therefore repre- 
sents an extremely interesting option 

Indeed, a TFET is a gated p — i — n diode where the 
gate voltage modulates the position of energy bands in 
the intrinsic channel in order to control interband tunnel- 
ing between the p-doped source and the n-doped drain. 
When the device is ON, current is dominated by inter- 
band tunneling, which instead is inhibited when the de- 
vice is OFF. Small energy gaps and small effective masses 
can allow to achieve large ON current, comparable to that 
of mundane semiconductor field-effect transistors. 

An impressive advantage of TFETs is the steep de- 
pendence of the current on the gate voltage, which 
can yield a Subthreshold Swing (SS) much smaller than 
60 mV/decade, the minimum value achievable at room 
temperature with FETs. The very small SS also allows 
to implement adequate switches with very small sup- 
ply voltages, and therefore to operate at extremely low 
power. 

Finally, carbon channels very often have the same dis- 
persion relations and mobility characteristics for elec- 
trons and holes, which is desirable for the optimization 
of complementary logic gates in terms of switching speed 
and power consumption 0]. 

TFETs based on carbon nanotubes [3, graphene 
nanoribbons 0, 0] and bilayer graphene J^TlgJ have been 
addressed theoretically, with simulations and analyti- 
cal models. Carbon nanotubes and graphene nanorib- 
bons are interesting in terms of the achievable perfor- 
mance, but device operation is extremely sensitive to the 
bandgap, which in turn strongly depends on the exact 
number of carbon atoms per ring or along the ribbon 
width, that is hardly controllable. Bilayer graphene is 
a two-dimensional channel which does not require very 
advanced lithography and provides a gap tunable with 
the vertical electric field, up to 0.15-0.2 eV. Satisfactory 



device operation is achieved with a supply voltage Vdd 
as small as 0.1 V [5j. 

Very interestingly, epitaxial graphene grown on the Si- 
terminated face of a SiC substrate is a two-dimensional 
material exhibiting a bandgap of 0.26 eV [7H9( and 
promising mobility [Toj j. Such semiconducting gap can 
be adequate to TFETs suitable for large-scale integrated 
circuits, but is not sufficient for FETs [llj . 

Here we consider the TFET illustrated in [Ha), in 
which the channel consists of a single graphene layer 
grown on a SiC substrate of overall thickness ts u b = 
100 nm. The top gate is separated by a silicon oxide 
layer t ox of 1 nm and is at voltage Vq. Following the 
evanescent mode approach [l2T ] , we can estimate the nat- 
ural variation length A of the electrostatic potential in 
the channel as A = ^(t ox + t c h), where t c h — 1.1 nm 
is the effective distance between the graphene layer and 
the ideal interface of the top dielectric layer [l3j]. Here, 
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FIG. 1: (a) Sketch of a TFET device based on graphene on 
SiC. (b) Band profile along the transport direction in the OFF 
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we have assumed the evanescent part of the potential to 
be dominated by the vicinity of the top gate, neglecting 
the backgate, and, second, we have assumed the same 
effective dielectric constant for both oxide and graphene 
regions. Because the gate length L c is always much larger 
than A, a long-channel behavior is expected. Therefore 
we write <fi(x) , the electrostatic potential on the graphene 
channel, in the form 

(x-L c )/\ | i 



q<p(x) = qA s e~ x/x + qA 



(1) 



where 4> c is the potential in the channel away from the 
contacts, q the absolute electron charge, As = E c (0) — 
E c (L c /2) and A D = E C (L C ) - E c (L c /2) respectively, 
where E c {x) represents the conduction band edge in the 
longitudinal direction x. 

For graphene on SiC we adopt the tight-binding Hamil- 
tonian proposed in [7[ , which accounts for the energy dis- 
persion curves 
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(2) 



of conduction (+) and valence (— ) bands, whose edges 
respectively give E c (x) and E v (x). Where m — 0.13 eV 
is an effective parameter providing an energy gap of 2m, 
t = 2.7 eV is the nearest-neighbor hopping energy for 
graphene and f(k) is the well known off-diagonal element 
of the 2x2 graphene p z -Hamiltonian. We assume the po- 
tential to be sufficiently smooth to rigidly shift graphene 
band edges. The band edge profiles for the device under 
investigation in the OFF state and in the ON state are 
shown in Fig. (Ub) and[IJc). 

The device current / is the sum of three major contri- 
butions: thermionic, interband tunneling and gap tun- 
neling (due to tunneling through the channel in the gap 
spectral region). Thermionic and gap tunneling currents 
degrade transistor performance, dominating the current 
when the device is in the OFF state, and therefore must 
be minimized. The channel current per unit width is 
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(3) 

where the integration runs over the Brillouin zone for con- 
fa: ) 

duction and valence bands E\ , calculated at the center 

k,± ' 

of the channel (at x — x c = L c /2). f is the Fermi-Dirac 
distribution function, v x is the component of the group 
velocity along the channel. The tunneling coefficient T 
is calculated via Wentzel-Kramers-Brilluoin approxima- 
tion in the tunneling regions at drain, at source or in the 
gap, as shown in Fig. [ljb). In tunneling regions, for a 
fixed energy E of the incoming particle, we solve the dis- 
persion curve ([2]) for the complex wavevector describing 
the evanescent mode. The integration of the evanescent 
mode path through the barrier gives T, as described in 
more detail in Ref. (Tlj . 

In order to compute the channel potential we need to 
self-consistently solve the vertical electrostatics: 

Q = -Cg{Vg — Ag — 4> c ) - C Su b(Vsub - A Sub - 4> c ) 
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FIG. 2: Transfer characteristics of a graphene-on-SiC TFET 
calculated for different parameters of the model. We analyzed 
the effects of a variation of doping concentration a in (a), of 
the supply voltage Vbs in (b), of the channel length L c in (c). 
In (d) we show an example as the TFET total current is due 
to the summing up of interband tunneling, in-gap tunneling 
and thermionic current contributions. 



Q = Qh + Qe 



(4) 



where C G = e x/t x, C Su b = esic/tsub are the top gate 
and the backgate capacitance per unit area, respectively, 
whereas Ag and A$ u b are the flat band voltages of the 
top and back gates (set here to V). The charge per 
unit area in the channel, due to holes and electrons {Qh 
and Q e ), can be also obtained by summing over their 
steady-state distribution in conduction band 
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where we included the possible tunneling barrier at 
source, drain or in-gap tunneling. 

Let us assume the same doping density for source and 
drain. We can therefore define the energy difference be- 
tween conduction(valence) edge and source(drain) Fermi 
energy at source(drain) S = E c (0) — (is = — — E V (L C ). 
It is clear from Fig. [ljb,c) that it is possible to shut 
off the device only if the drain-to-source voltage satisfies 
I Vbs | < (Egap — 25) I q, which in the limit of weak doping 
becomes \V DS \ < 0.26 V. 

Ideally, when the previous condition is satisfied, mini- 
mum device current is obtained when the channel is in- 
trinsic at L c /2, and therefore has zero net charge. Ac- 
cordingly the channel mid gap is — q<f) c = (ps + /ir>)/2 
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FIG. 3: Jon and /off table for TFET devices with Vdd = 
0.1, 0.15, 0.2, 0.25 V, and L c between 10 and 50 nm. We 
consider here t ox — 1 nm and a — 3 x 10 16 m -2 . A dashed 
line marks the field where /on/^off > 10 4 . 



and from ^ we find the OFF-state gate voltage Voff = 
Vds/2. The OFF current /off is defined as the cur- 
rent obtained with Vqs = Voff arid Vds = Vdd- The 
ON current Jon is the current obtained with Vgs = 
Voff + V D d, and Vds = Vdd- 

In Fig. [2] we show the TFET transfer characteristics 
calculated for different values of the main parameters of 
the model: the doping fraction a (which is a monotonous 
function of 6), Vds, the channel length L c . In Fig. [2)(a) 
we show the effect of doping concentration a = 10 16 , 
10 17 , 2 x 10 17 , 3 x 10 17 m~ 2 , corresponding respectively 
to 6 « 5, 40, 65, 90 meV. 

As expected, the increase of 5 leads to a thinning of 
the range of the OFF-state plateau, in which current is 
essentially reduced to the thermionic contribution. An 
initial increase of S pushes down the source conduction 
band edge with respect to (Xg, and up the drain valence 
band with respect to fin, reducing the thermionic hole 
and electron currents respectively. For larger a (and 5), 
one cannot fully shut the channel off because |Vos| + 2(5/g 
becomes larger than E gap . 

In Fig. [3](b) we analyze the operation of the TFET 
for different values of the source-drain voltage Vds ■ The 
OFF-state plateau is reduced with increasing Vds, /off 
is constant because limited to the thermionic contribu- 
tion, whereas larger /on is obtained. For Vds ~ 0.25 V 
and above it is not possible to shut the TFET off any- 
more. 

In view of transistor scaling we analyze the effect of a 
variation of the channel length L c (Fig [21(c)). We show 
that for L c > 30 nm the OFF current does not depend 
on L c and is determined by the thermionic current. How- 
ever, a further reduction of L c leads to an increase /ofFj 
due to an increased transparency of the in-gap tunneling 
barrier. The oxide thickness affects tunneling through 
the contacts (A oc t ox ), with the main effect of exponen- 
tially reducing the tunneling log/ 1 oc — t ox . In Fig. [2](d) 
we show an example of the three current contributions 
occurring in the device. In ON state the tunneling cur- 



rent dominates the device, while it is suddenly suppressed 
going towards the OFF state, where thermionic current 
contribuites to the subthreshold swing and evantually in- 
gap tunneling determines the OFF state current plateau. 

An important performance measure of a transistor is 
the obtainable /on //off ratio, therefore in Fig. we 
plot /on against /off- As noted also before, to shut the 
device off and have a sufficent /off dynamics it is re- 
quired Vds < 0.2 V, togheter with sufficiently long chan- 
nel (L c > 30 nm) in order to limit in-gap tunneling. On 
the other hand the /on increases with Vdd- Therefore 
the best performance is obtained for Vds = Vdd = 0.2 V. 
Indeed we note that an /on//off ratio exceeding 10 4 , 
considered adequate for low power operation 14| , can be 
obtained with a supply voltage of 0.15 and 0.2 V, and 
with a channel length L c not smaller than 30 nm. 

Now we analyze the performance of our device for dig- 
ital applications. In Fig. 0Ja) we plot the transfer char- 
acteristics of a CMOS inverter with identical n and p- 
type TFETs, with the channel of graphene on SiC with 
L c = 30 nm and t ox = 1 nm, for a supply voltage of 0.1 V 
and 0.2 V. We extract the noise margins (NMs) for the 
CMOS inverter, which are a measure of logic gate robust- 
ness to disturbs and noise, as NM = 22 mV and 37 mV, 
respectively for Vdd =0.1 and 0.2 V. NMs are not ex- 
tremely good, but in line with what one could expect, 
given the low supply voltage. 
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FIG. 4: (a) Voltage transfer curve for a CMOS inverter based 
on the present device and operating at Vds = 0.1 and 0.2 V. 
Dashed lines highlight the points of the curve in which the 
slope is equal to unity, (b) Cutoff frequency calculated for 
TFETs with Vds = 0.1 and 0.2 V, for L c = 30 and 50 nm. 

We also performed an analysis of transistor intrin- 
sic delay and transition frequency fx of the considered 
TFET. We consider operating voltages of 0.1 V and 0.2 V 
and a channel length of 30 and 50 nm. In particular in 
Fig. |U(b) we show the cutoff frequency as a function of 
the overdrive, calculated in the quasi-static approxima- 
tion [H /t = g m /{2irC G L c ), where g m = dI D s/dV GS is 
the transconductance, and Cq = dQ /dVos evaluated for 
Vgs — Vdd + Voff, is the capacitance seen by the gate 
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contact. The intrinsic delay time, which represents a fig- 
ure of merit for digital application, is t = VddCg/^ON- 
For a supply voltage of 0.1 V, we obtain r = 0.33 ps 
and 0.54 ps, respectively for a channel of length 30 and 
50 nm. With a supply voltage of 0.2 V, r = 0.12 ps for 
L c = 30 nm and r = 0.22 ps for L c = 50 nm. 

In conclusion, we examined the maximum achievable 
performance of a TFET based on epitaxial graphene on 
SiC substrate, by means of an accurate semi- analytical 
model. We explored the parameter space analyzing its 
performance under variations of the doping level, supply 



voltage, channel length and oxide thickness. We show 
that an Ion/^off exceeding 10 can be easily achieved 
by choosing adequate geometry and bias, and a channel 
length not smaller than 30 nm, even for a supply volt- 
age as small as 0.15 V. Noise Margins of a minimum size 
inverter are reasonable but not impressive. Intrinsic de- 
lay times and transition frequency are very promising, 
especially considering the very low supply voltage, and 
make the proposed TFET a strong candidate for extreme 
low-power carbon nanoelectronics. 
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